function A = symm(A)
% A = symm(A)
% Returns the symmetric part of A

  A = 0.5 * (A + A');

end